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Ballooning instabilities are investigated in three-dimens- 
ional magnetic toroidal plasma confinement systems with low 
global magnetic shear. The lack of any continuous symmetry 
in the plasma equilibrium can lead to these modes being lo- 
calized along the field lines by a process similar to Anderson 
localization. This produces a multibranched local eigenvalue 
dependence, where each branch corresponds to a different unit 
cell of the extended covering space in which the eigenfunction 
peak resides. These phenomena are illustrated numerically 
for the three-field-period heliac H-1, and contrasted with an 
axisymmetric s-a tokamak model. The localization allows a 
perturbative expansion about zero shear, enabling the effects 
of shear to be investigated. 

PACS numbers: 52.35. P, 52.55.H, 71.23.A 



Ballooning instabilities are pressure-driven ideal mag- 
netohydrodynamic (MHD) instabilities which limit the 
maximum (3 (plasma pressure/magnetic pressure) that 
can be obtained in a plasma. They are localized about 
regions where the field lines are concave to the plasma, 
which are known as unfavourable regions of curvature. 
Another localizing influence is the magnetic shear, which 
measures the rate at which neighboring field lines at dif- 
ferent minor radii separate as they wind their way around 
the torus. Large shear helps stabilize these modes, 
thereby playing an important role in the MHD stabil- 
ity. In this paper however we consider the effects of very 
small or zero shear, such as occurs in the heliac class of 
stellarators or in the shear-reversal layers of an advanced 
tokamak. 

We begin by making the usual assumption that the 
magnetic field lines map out nested flux surfaces, or mag- 
netic surfaces. These are labeled using a normalized- 
toroidal-flux variable s, which varies between zero at 
the center of the plasma and unity at the plasma edge. 
Within each surface the straight-field-line poloidal 9 and 
toroidal £ angle variables are defined such that the field 
lines appear as straight lines in the (9, £) plane. The 
magnetic field may then be written B = V^xV?/; — 
qVOxVip = VaxV-0, where the field-line label a = 
C — q9. Here, 2Tvt/j represents the poloidal magnetic flux, 
while q — q(s) is the safety factor (inverse of rotational 
transform), which is equal to the average number of 
toroidal circuits traversed by a field line per poloidal cir- 



cuit traversed around the torus. 

Ballooning modes can be characterized as having a 
long parallel and short perpendicular wavelength with 
respect to the field lines. By ordering the perpendicular 
wavelength to be small and expanding to lowest order in 
an asymptotic series the local mode behavior can be ex- 
pressed by a one-dimensional equation along a field line 
jjj. Taking the plasma to be incompressible, the balloon- 
ing equation may be written 



(1) 



where the eigenfunction £ is related to the mode dis- 
placement while the eigenvalue A is equal to the mode 
growth rate squared. This represents the local stability, 
local to a field line. In forming global modes, ray trac- 
ing must be performed in the three-dimensional A phase 
space to determine which of these local solutions also sat- 
isfies the poloidal and toroidal symmetry requirements on 
the phase of the mode || . This paper is concerned with 
the functional dependence of A on its three arguments. 

The ballooning coefficients are functions of local pa- 
rameters and (9 — 9k) "secular" terms, where the secular 
terms represent the effects of magnetic shear 
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where B represents the field strength, p is the 
plasma pressure, and the Jacobian is given by J = 
(VVj'V^X VC) • The local integrated shear is given 
by TZ + 9 d^q = — Va-Vi/>/| V-0| 2 , and the normal and 
geodesic components of the magnetic curvature vector 
k = e||-Ve|| (where e|| = B/B) are given by n n = 
k-Vs/\Vs\ and n g = K-VsxB/\BVs\, respectively [Q. 
The parameter 9k is related to the direction of the mode 
wave vector. 

The coefficient A is positive definite, allowing the bal- 
looning equation to be transformed into the Schrodinger- 
like form [d 2 /d6» 2 +E- V] A 1/2 £ = 0, where the "po- 
tential" 
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and E = —\J 2 so that instability occurs whenever the 
"energy" is negative ||. 

We first consider the case of a two-dimensional equi- 
librium, which contains an ignorable third dimension. 
Specifically we use the s-a tokamak model, which de- 
scribes an axisymmetric equilibrium analytically by as- 
suming circular flux surfaces and a large aspect ratio || . 
Being axisymmetric, the toroidal angle is ignorable while 
the local parameters are 2-7r-periodic in the poloidal an- 
gle. The potential may be written 
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where 5T = d(ln q) /<9(ln r) is a measure of the global shear 
across the magnetic surface, a = — 2Rq 2 d r p/ B 2 is a mea- 
sure of the pressure gradient, h = s"(0 — 6k) — a sin rep- 
resents the integral of the local magnetic shear along the 
field line, and hi = Ah/ AO is equal to the local shear (we 
use s~ and a instead of the usual labels s and a in order 
to avoid confusion with the normalized-toroidal-flux vari- 
able s and the field-line label a). Here, r and R represent 
the minor and major radii of the torus, respectively. The 
energy is given by E = —XRq 2 /B 2 . 
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FIG. 1. Plots of V(6) - E vs. 6 (solid lines), for the s-a 
model at 9k = and the most unstable a for s" — 0.1, 1. 
Dashed lines represent the wavefunctions _4 1//2 £. 

Two examples of V(0) — E for the s-a model are 
shown in Fig. 0, along with the corresponding "wave- 
functions" A 1 / 2 ?;, which are very similar to the eigenfunc- 
tions £ for these cases. These correspond to a low-shear 
(s = 0.1) and a high-shear (s = 1) case. The wavefunc- 
tion (and eigenfunction) peaks occur at multiples of 2n 
in 9, where the normal curvature component (given by 
K n = — R^ 1 cos 8) is most unfavourable. This behavior is 
modulated by the secular terms, which are proportional 
to s" and have the effect of localizing the eigenfunction 
around 9 « 9k- 



The first term in Eq. (||) is alternatively stabilizing and 
destabilizing as varies, being proportional to both the 
normal curvature and the pressure parameter. This term 
is modulated by a function of the integrated shear. In the 
large-|6*| limit the influence of this shear causes V — > 0, 
ensuring that all unstable solutions (both the wavefunc- 
tion and eigenfunction) are exponentially localized. Con- 
sequently, the s-a model is stable to the more extended 
interchange modes. The second term is purely stabilizing 
and represents the effects of the local shear. For large \9\ 
this term is 0(9~ 4 ) compared with 0(9~ 2 ) for the first 
term, so the effects of shear are actually more dominant 
in the first term, with this term playing a significant role 
at more moderate values of 9. 
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FIG. 2. Eigenvalues A vs. Ok, for the same s", a as Fig. 
The poloidal branch labels n are marked for the second case. 

In Fig. H we investigate the (^--dependence of the local 
eigenvalues. This is periodic after a distance 2-7T in as 
expected from the invariance of the ballooning equation 
under the poloidal mapping operation, P : 9 h=- 9 + 2-7T, 
Ok l— * Ok + 2tt, a t— > a — 2nq. An increase in Ok by 27r is 
therefore associated with an eigenfunction shift of 2ir in 
the poloidal angle along the field lines. This enables us to 
define the poloidal branch label n, representing the 2-7T in- 
terval in where the eigenfunction is at a maximum. Val- 
ues of this branch label are marked in Fig. ||(b), where n 
is defined such that the eigenfunction peak occurs around 
w 2im. 

The high-shear case gives well defined branches, in con- 
trast to the low-shear case where the branches merge so 
that there is no stable interval of Ok ■ The difference is due 
to the mode structure along the field lines (see Fig. |l|). 
As the shear is decreased the eigenfunctions become more 
extended along the field lines and the difference between 
eigenfunctions of neighboring branches is decreased. In 
the s — > limit all branches merge and the most un- 
stable eigenfunction becomes periodic after a distance 
2tt in 0, due to the vanishing of all secular terms. The 
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generalized solutions of the ballooning equation are then 
Block waves, which can be written in the Floquet form 
£ = <pk{9) exp(iK0), where K represents the Block wave 
vector and (pa is a 27r-periodic function in 9. 

In contrast to the above two-dimensional case, a gen- 
eral three-dimensional equilibrium contains local param- 
eters which exhibit only quasiperiodic variation along the 
field lines, since, when q is irrational, a field line will 
cover a magnetic surface ergodically. We argue that this 
breaking of continuous symmetry leads to the ballooning 
eigenfunctions being localized along the field lines, even 
in the absence of magnetic shear. 

We study an equilibrium which has the standard mag- 
netic configuration of the three-field-period H-l heliac 
at the Australian National University Q. The pressure 
profile is similar to that which produces marginal stabil- 
ity throughout the plasma ||, but increased by 20% so 
that the plasma becomes highly unstable with a /3 of 1% 
instead of the marginally stable value of 0.8%. The pre- 
conditioned VMEC code H was used to calculate the equi- 
librium and then a mapping code was used to transform 
it into straight-field-line Boozer coordinates. In doing 
this the poloidal and toroidal angles were defined such 
that (9, () = (0, 0) corresponds to a symmetry point on 
the outer side of the stellarator about which stellarator 
symmetry holds [ jTo| . This point occurs deep within a 
region of unfavourable curvature. 
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FIG. 3. Safety factor q vs. surface label s, for the H-l con- 
figuration studied. 

A plot of the safety factor variation is shown as Fig. [|. 
The shear is small throughout the plasma volume with 
q varying by less than 5%. The shear parameter of the 
s-a model is approximately given by s~ « 2sd s q/q. Us- 
ing this expression, the maximum absolute value of s for 
this configuration is only around 0.12, and occurs at the 
plasma boundary. 
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FIG. 4. Contours of local eigenvalues for the H-l configu- 
ration at 6k = 0. Areas labeled by 's' represent stable regions. 
Some corresponding eigenfunctions are shown in Fig. 0. 
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FIG. 5. Eigenfunctions £ vs. poloidal 9 and toroidal £ an- 
gles, at the points marked by 'x' (two eigenfunctions), '□' 
and 'A' in Fig. ^. The branch labels (m,n) are marked. 



A representation of the phase-space dependence of the 
(most unstable) local eigenvalue for this configuration 
with 6k = is shown as Fig. ^, with some corresponding 
eigenfunctions in Fig. ||. Despite the fact that the shear is 
small (s w —0.06 at s — 0.70), the eigenfunctions are well 
localized along the field lines and are mostly confined to 
a single 2ir x 2ir/M unit cell of the (9, £) covering space, 
where M represents the number of identical field periods 
in the stellarator (M = 3 for H-l). We use this to de- 
fine the poloidal and toroidal branch labels, which identify 
the unit cell of the covering space where the eigenfunc- 
tion has a maximum. Each of the cases shown contains 
well defined branch labels, due to the strong localization 
of the eigenfunctions. This is in sharp contrast to the 
more extended eigenfunctions of the s-a model at a simi- 
lar value of |s| (see Fig. 0(a)), suggesting that something 
other than shear is responsible for the localization. This 
is confirmed by setting d^q = in Eq. (0) to eliminate 
the secular terms altogether. Once again, we find the 
eigenfunctions to be well localised, and indeed to be al- 
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most identical to their counterparts in Fig. M. 
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FIG. 6. Same as Fig. [l], except that this plot corresponds 
to the (a) '□' and (b) 'A' eigenfunctions of Fig. ^. The 
small-scale fluctuations in V(6) are a consequence of ripple 
effects from the 36 toroidal held coils. 

This localization can be attributed to Anderson local- 
ization, the process by which electron wavefunctions be- 
come localized in space due to the presence of impurities 
in an otherwise perfect crystal In this case how- 

ever it is the quasiperiodicity of local parameters, due to 
the incommensurate periods of the toroidal and poloidal 
variations on a held line when q is irrational, which cause 



the localization [fLz 



13. We plot two examples of the 



"potential" in Fig. |€ , which is quasiperiodic in the zero- 
shear limit. From this plot it is not immediately clear 
where the wavefunction A 1 / 2 !; will reside, since potential 
"wells" exist at many places along the field lines. 

The wave packet location and the resulting eigenvalue 
dependence can be understood in terms of the normal 
curvature. First consider the (m, n) = (0, 0) branch, 
labeled by '□' in Figs. ^ and p§. This eigenfunction 
contains a maximum around the {9, £) ~ (0, 0) region 
of unfavourable curvature, corresponding to a w in 
Fig. [l|, independent of s. The (m, n) = (3, 1) branch on 
the other hand appears to be most unstable around the 
a ps 2%(m/M - nq) ps 0.76 field line when s = 0.70, 
this being the field line that passes through the same 
(6, C) ~ (0, 0) unfavourable region, but only after travers- 
ing the torus once in both the poloidal and toroidal di- 
rections. Similarly, other branches contain eigenfunc- 
tions which peak around the same region of unfavourable 
curvature, but only after the field line has undergone 
m/M toroidal and n poloidal transits around the torus. 
The most-unstable field line of each branch with n ^ 
will therefore be a function of q(s), resulting in a com- 
plex "multibranched" eigenvalue structure, such as that 
shown in Fig. ^. 

While this behavior is qualitatively correct, in practise 
we find a small discrepancy between the predicted and 



observed positions of the maximum eigenvalue of each 
branch juj. For example, the (3,1) branch is actually 
most unstable at a w 0.85 for s = 0.70, in comparison 
with the predicted value of a ps 0.76. This "shifting" 
effect can be investigated using a two-dimensional per- 
turbative expansion of the ballooning equation, where we 
first solve along the most-unstable field line ao with the 
global shear set to zero, and we then treat finite (a — ao), 
d^q as a perturbation. 

Provided q is sufficiently irrational, the zeroth-order 
eigenfunction estimate will be square integrable and will 
be mostly confined to a single unit cell of the covering 
space. We perturb this solution by including terms in 
both d^q and (a — ao), which are assumed to be of the 
same order, and we consider all branches by including the 
branch labels explicitly. The final eigenvalue estimate to 
second order can then be written in the compact form 
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where e = d^q, the A^ coefficients are functions of the 
magnetic surface only, and for this H-l configuration 
ao = so all functions are of definite parity. Each of 
the above terms has a simple physical interpretation. 
The first represents the unperturbed estimate, which is 
corrected by the following two terms which act on all 
branches equally. These are stabilizing and represent 
the effects of shear on the (0, 0) branch. The fourth 
term is also stabilizing and represents the effects of shear 
as a field line is followed through one or more com- 
plete poloidal rotations. Finally there is the last term, 
which contains the field-line eigenvalue dependence and 
the eigenvalue peak "shift" term. This shift can be in- 
vestigated by ordering the eigenfunction to be well lo- 
calized along the field line, so that all odd-parity terms 
vanish at leading order. To lowest order we obtain 
Tl + {d^q){0 — 9k) ps at the mode peak, showing that 
the most unstable part of the branch occurs where the 
integral of the magnetic shear approximately vanishes at 
the mode peak. 

We have implemented this expansion numerically and 
have found that in most cases the expansion provides 
good agreement with eigenvalues calculated directly, ex- 
cept near the plasma edge and around the q = 8/9 ratio- 
nal surface (containing 67r-periodic local coefficients). In 
particular, the predicted "shifts" in the position of the 
eigenvalue maximum are in approximate agreement with 
those observed directly and those obtained by assuming 
that 1Z + {d 1 pq){9 — 9k) ~ at the mode peak, support- 
ing our interpretation of the mechanism behind this term 
(details of these results will be reported elsewhere). 

The numerical calculations were performed on the Aus- 
tralian National University Supercomputer Facility's Fu- 
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jitsu VPP300 vector processor. We wish to thank Dr. 
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